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Abstract 

Self-organization is pervasive in development, from symmetry breaking in the early embryo 
to tissue patterning and morphogenesis. For a few model systems, the underlying molecular 
and cellular processes are now sufficiently characterized that mathematical models can be 
confronted with experiments, to explore the dynamics of pattern formation. Here, we 
review selected systems, ranging from cyanobacteria to mammals, where different forms of 
cell-cell communication, acting alone or together with positional cues, drive the patterning 
of cell fates, highlighting the insights that even very simple models can provide, as well as 
the challenges on the path to a predictive understanding of development. 


Introduction 

During development, an elaborate spatial structure is gradually established, through the 
allocation and spatial rearrangement of different cell types. Although much of development 
can be understood as a cascade of successive refinement steps, with early-formed 
landmarks providing positional cues for further subdivisions (Wolpert, 1969), some of these 
steps occur in the absence of overt positional cues, as in the development of the early 
mouse embryo (White et al., 2018). Self-organization is also evident in the ability of 
previously dissociated cells to recapitulate essential steps of development in vitro, from 
germ layer and axial patterning in cultures of mammalian stem cells (Beccari et al., 2018; 
Warmflash et al., 2014), to the development of organoids from cell aggregates (McCauley 
and Wells, 2017; D. A. Turner et al., 2016), to traveling waves of gene expression in cultures 
of presomitic mesoderm cells (Tsiairis and Aulehla, 2016). Self-organization, which may be 
defined broadly in the context of development as the spontaneous emergence of spatial 
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order through cell-cell interactions, mediated e.g. by diffusing ligands or mechanical forces, 
can take different forms (Sasai, 2013). Here, we focus on self-organized tissue patterning, 
understood as the spontaneous patterning of cell fates. In this sense, self-organization is to 
be distinguished from self-assembly, as in the formation of epithelial structures such as 
tubules (Montesano et al., 1991) or somites (Dias et al., 2014), or the sorting of cells that 
have already acquired different fates through differences in cell adhesion or 

contractility (Fagotto, 2014); and from self-organized morphogenesis, which is not 
necessarily linked to fate patterning, as in the development of branched organs (Lu and 
Werb, 2008; Ochoa-Espinosa and Affolter, 2012), and in tissue buckling (Shyer et al., 2013). 
And while our focus is on the scale of the tissue, development can also involve self- 
organization within individual cells (Halatek et al., 2018), as in the polarization of the early C. 
elegans embryo (Goehring et al., 2011). 


Principles and models of self-organized patterning 

The dynamics brought about by mutual cell-cell interactions are difficult to anticipate, and 
the experimental study of self-organization is inseparable from the exploration of 
mathematical models for pattern formation. Indeed, while some form of self-organization 
was implicit in the earlier concept of embryonic regulation (Roth, 2011), an explicit 
mechanism for spontaneous symmetry breaking - for the spontaneous generation of form - 
during development first came with the reaction-diffusion models introduced by 

Turing (Turing, 1952), in which cells interact through diffusing chemicals (Box 1). In a simple 
reaction-diffusion model, an activator stimulates its own production and the production of 
an inhibitor; if the inhibitor diffuses faster than the activator, a stable pattern can emerge 
from a uniform background through the amplification of small differences. In a two- 
dimensional geometry, which is a commonly adopted to reflect the geometry of tissues such 
as epithelia or skin, different choices of model parameters give rise to regular arrays of 
stripes or spots, which have notably been likened to animal coat patterns (Koch and 
Meinhardt, 1994; Murray, 2003). Beside stable patterns, reaction-diffusion models can also 
give rise to oscillations and traveling waves, and have been proposed to describe 
somitogenesis (Cotterell et al., 2015). In addition to a large body of theoretical work (Koch 
and Meinhardt, 1994; Murray, 2003), Turing's ideas have stimulated the search for signaling 
molecules that could mediate pattern formation, and candidate activator-inhibitor pairs 
have been proposed to underlie patterning in diverse model systems, from hair and feather 
primordia (Michon et al., 2008; Sick et al., 2006) to palate rugae (Economou et al., 2012). 
Yet, as popular as the activator-inhibitor paradigm remains, models that are more closely 
tied to mechanism often require elaborations on this simple scenario, like the inclusion of 
non-diffusible species such as transcription factors, encoding a cell state that responds to 
signals and governs their production. And as Turing himself recognized, diffusing chemicals 
are just one way that cells can communicate. Long-range signaling can also be mediated by 
membrane-bound ligands carried by cellular projections, or by mechanical forces. 
Additionally, motile cells can arrange themselves into patterns of cell density, which may 


couple with cell fate. As mathematical models have explored, spontaneous cell aggregation 
can occur through secretion of an attractant (Keller and Segel, 1970), or through mechanical 
cues transmitted by the extracellular matrix (Murray et al., 1983). Although, strictly 
speaking, Turing models involve interactions based on diffusing chemicals, it is noteworthy 
that models assuming different forms of cell-cell communication can be very similar in their 
mathematical structure and patterning behavior. And to some extent, different interaction 
modalities can be subsumed under an generic "interaction function" (Hiscock and Megason, 
2015) (Box 2), which could provide a common framework to address the logic of pattern- 
forming systems. For instance, pattern-forming models commonly involve local activation 
and long-range inhibition, but, as we discuss below, whether this must always be the case 
remains a matter of debate. 


In addition to involving different forms of interactions, self-organized patterning can give 
rise to different forms of spatial structure. As a general rule, the range of cell-cell 
interactions, e.g. the diffusion range of signaling molecules, governs the size of pattern 
elements, determining different classes of patterns and calling for different modeling 
approaches. Long-ranged interactions give rise to smoothly varying patterns of signaling 
activity and gene expression, at least before a stage where fate boundaries sharpen, and the 
relevant variables are commonly described as continuously varying functions of space. 
Short-ranged interactions, on the other hand, can form patterns of isolated cells, like the 
salt-and-pepper arrangements of neural cells arising through lateral inhibition, and call for 
models that explicitly describe individual cells. Depending on how the range of interactions 
compares to the size of the tissue, the output of self-organization can range from tissue- 
wide gradients, defining a positional cue for downstream patterning events, as in left-right 
patterning in the mouse (Nakamura et al., 2006), to patterns comprised of a small, invariant 
number of elements, such as digits in mammals (Raspopovic et al., 2014) and stripes in fishes 
(Nakamasu et al., 2009), to a uniform covering with regularly spaced elements, such as hair 
and feather primordia (Michon et al., 2008; Sick et al., 2006). 


Self-organized patterning is commonly contrasted with the specification of cell fates by 
positional cues, yet it is likely that much of development involves an interplay between pre- 
existing cues and cell-cell interactions, more than one or the other (Green and Sharpe, 
2015). Purely self-organized patterns, arising from initial or temporal fluctuations, are by 
essence variable, and models for disordered patterns such as hair primordia can compare 
with experiment only in an average sense, judging e.g. by the typical size or spacing of 
pattern elements in control or perturbed conditions. But where the output is a stereotyped 
arrangement like digits in mammals (Raspopovic et al., 2014), or a gradient with an invariant 
orientation, like Nodal in zebrafish (Muller et al., 2012), it must be that particular initial or 
boundary conditions steer patterning towards a reproducible outcome. Such a bias may 
arise from pre-existing positional cues, or from the geometry of the tissue itself, possibly 
including its growth. Reaction-diffusion models that would generically produce variable 


outcomes follow a determined sequence in a growing domain (Crampin et al., 1999), and the 
borders of stem cell colonies grown on patterned substrates, acting as sinks for diffusing 
ligands, have been proposed to provide a positional cue for stereotyped radial 

patterning (Etoc et al., 2016). Such a role for initial or boundary conditions is not restricted 
to patterns of a few elements, since large-scale patterns that are constrained to form 
sequentially can exhibit long-range order, like the arrangement of ommatidia in the 
Drosophila eye, which forms in the wake of a traveling differentiation front (Lubensky et al., 
2011; Roignant and Treisman, 2009), or the arrangement of feathers in chick (Ho et al., 
2019). And in the same way that self-organized patterns never emerge on a perfectly blank 
slate, there is much evidence to suggest that cell-cell interactions can refine the 
interpretation of morphogen gradients, or even that the receiving cells contribute to shaping 
the gradient, as in the patterning of the vertebrate neural tube by Sonic Hedgehog 

(Shh) (Ribes and Briscoe, 2009). The establishment of gradients of Nodal and its inhibitor 
Lefty in the zebrafish embryo is one rare example where molecular diffusivities were 
quantified and supported a Turing model (Muller et al., 2012), yet recent experiments 
suggest that inhibitory feedback from Lefty is not strictly required for patterning, since 
uniform Lefty expression rescued lefty mutants (Rogers et al., 2017) - thus many gradient- 
forming systems may exist on a continuum between positional information and self- 
organized patterning. 


As this overview should make clear, self-organized patterning encompasses very diverse 
phenomena and model systems. Rather than aim for a comprehensive survey, we have 
chosen to focus primarily on patterns of repeated elements, and on a few experimental 
settings in which experimental data can now be confronted with models in a semi- 
quantitative manner, ranging from a simple one-dimensional system in prokaryotes to the 
more complex two-dimensional pattern of the lizard skin. For reasons of space, many other 
systems of interest are not covered here, including the establishment of tissue-wide 
gradients in vertebrates through Nodal and Lefty (Muller et al., 2012; Nakamura et al., 2006) 
and during planarian regeneration (Stückemann et al., 2017); somitogenesis (Hubaud et al., 
2017; Pourquié, 2011; Sonnen et al., 2018); and patterning in plants (Bhatia and Heisler, 
2018). Self-organized patterning in the early mammalian embryo has been covered in a 
separate review (White et al., 2018). The systems covered here illustrate various forms of 
cell-cell communication, from diffusing ligands to mechanical forces, different modeling 
approaches, from continuous to cell-based, and various forms of interplay between 
positional cues and self-organization. As these examples also illustrate, and as we further 
discuss in closing, self-organized patterning often involves an interplay between multiple 
signals, whether chemical or mechanical, of which a subset may suffice to explain the 
observed outcome. While such an apparent redundancy can be rationalized as a conferring 
robustness on development, it implies that striking a balance between model complexity, 
predictive power, and interpretability is a delicate exercise, and that proving or disproving 


any particular mechanism may require stringent tests, beyond assessing the effect of 
uniform perturbations on the final pattern. 


A one-dimensional self-organizing system in cyanobacteria 

A comparatively simple instance of a self-organized pattern is observed in several species of 
cyanobacteria such as Anabaena, which form colonies consisting of one-dimensional 
filaments with a periodic pattern of two distinct cell types (Kumar et al., 2010; Meeks and 
Elhai, 2002; Wilcox et al., 1973a; 1973b). In the presence of nitrogen sources, the filaments 
comprise only mitotically active vegetative cells that carry photosynthesis. In contrast, under 
nitrogen-limiting conditions, some cells stop dividing and differentiate into heterocysts that 
can fix atmospheric nitrogen. A key feature of this system is that heterocysts appear at 
regular intervals, forming a one-dimensional pattern of heterocysts separated by about ten 
vegetative cells. This pattern is dynamic: as vegetative cells continue to divide, new 
heterocysts differentiate in the middle of the intervals between pre-existing heterocysts, so 
that a regular pattern is maintained over time. 


This basic phenomenology suggested that lateral inhibition from differentiated heterocysts 
maintains a regular spacing pattern. An early mathematical model that postulated a diffusing 
inhibitor, with an exponential decay away from heterocysts, agreed well with the measured 
distribution of intervals between the first-formed heterocysts (Wilcox et al., 1973a; Wolk 
and Quine, 1975). Subsequent experimental studies provided support for this model. 
Heterocyst differentiation is positively regulated by a transcription factor, HetR, that is 
inhibited by two short diffusible peptides, HetN and Pats, that inhibit HetR (Corrales- 
Guerrero et al., 2013; Kumar et al., 2010; Risser and Callahan, 2009). While HetR and PatS 
display basal fluctuating levels under nitrogen-rich conditions (HetN only appears later in 
differentiated heterocysts and acts mostly in pattern maintenance)(Corrales-Guerrero et al., 
2015), an upstream regulator of HetR is activated upon nitrogen stress and elevated levels of 
HetR in turn up-regulate PatS expression (Figure 1). 


The identification of these molecular players, and the resulting mathematical models (Brown 
and Rutenberg, 2014; Mufioz-Garcia and Ares, 2016), have provided a clearer picture of 
patterning in cyanobacteria. Yet they leave some intriguing features of heterocyst 
differentiation unexplained. Early observations suggested that differentiating heterocysts do 
not emerge in isolation, as could be expected from a simple lateral inhibition mechanism. 
Rather, strings of differentiating cells are often observed, that later resolved into a single 
heterocyst (Wilcox et al., 1973a) (see also the pattern of HetR expression in Figure 1A). This 
suggested a "two-stage" model, where differentiation is first initiated in groups of 
competent cells, followed by their resolution (Corrales-Guerrero et al., 2015; Meeks and 
Elhai, 2002). The emergence of the groups at regular intervals may suggest spatially 
correlated fluctuations in the basal expression of HetR. Indeed, live imaging of a HetR-GFP 


fusion protein and analysis of spatial fluctuations under nitrogen-rich conditions suggest 
coupling between neighboring cells (Corrales-Guerrero et al., 2015). Alternatively, spatially 
correlated cell cycle phases, on a scale commensurate with the heterocyst pattern (Meeks 
and Elhai, 2002), and asymmetries arising through division in protein content (Risser et al., 
2012) or cell size (Mitchison and Wilcox, 1972), might modulate cell competence. When 
most mathematical studies of pattern formation emphasize the selective amplification of 
random fluctuations as the main driver of pattern selection, the example of Anabaena thus 
points to the possible role of preexisting structure within a population of equivalent cells. 
Such a pre-existing bias, resulting from a difference in birth order, has been proposed to 
underlie the AC/VU decision in C. elegans (Karp and Greenwald, 2003). In cyanobacteria, 
structured fluctuations on scales large and small could facilitate symmetry breaking in the 
two stages of patterning. Live imaging along with mathematical models that explicitly 
integrate correlated fluctuations should reveal how pattering builds on a structured 
background. 


Also, can one rationalize the existence of two diffusing inhibitors, HetN and PatS, where one 
might seem enough to support the emergence of the pattern and its maintenance during 
growth? Could it be that heterocyst selection among competent cells calls for a fast inhibitor 
turnover, while an inhibitor with a slower turnover can support lateral inhibition from 
differentiated heterocysts at a reduced "cost"? Or does the transition from mutual to lateral 
inhibition require a subtle regulation of inhibitor levels, that is more "readily" achieved by 
combining multiple inhibitors, with individually simpler regulation? By formalizing such 
questions, modeling could help define further experiments to dissect the regulatory 
architecture underlying the formation and maintenance of a regular spacing pattern in 
cyanobacteria. 


Stripes and dots: self-organization in flies 

While patterning in cyanobacteria affords a simple example of a one-dimensional spacing 
pattern, the suggestion that such arrangements arise from lateral inhibition first arose from 
the study of sensory organ patterns on the epidermis of insects, and Wigglesworth's 
observations on the development of sensory bristles in the blood-sucking bug Rhodnius. 
Rhodnius grows in size through five nymphal stages, with new sensory bristles emerging at 
each molt. Noting that these new bristles develop in the largest intervals between existing 
ones, Wigglesworth proposed that bristles exert an inhibitory influence on their 
surroundings, possibly by consuming a uniformly produced bristle-promoting factor 
(Wigglesworth, 1940), as in a substrate-depletion model (see Box 1). A regular spacing 
pattern would thus be maintained as the body enlarges. Four decades later, genetic and 
molecular studies in Drosophila led to revisit this model. In contrast with Rhodnius, the adult 
fly is produced at once to its final size by the differentiation of imaginal tissues through a 
complete metamorphosis, thus bristles develop in their definitive pattern, without 


subsequent insertion of new bristles. Where bristles can form is governed by the expression 
of proneural transcription factors such as Achaete and Scute, which endow cells with the 
competence to become sensory organ precursor cells (SOPs). SOPs are then selected from 
among the competent cells through lateral inhibition mediated by the Delta ligand and the 
Notch receptor (Ghysen and Dambly-Chaudière, 1989; Simpson, 1990), which antagonizes 
the activity and expression of the proneural factors through its E(spl)-HLH targets (Delidakis 
et al., 2014). The pattern of proneural gene expression thus largely prefigures the layout of 
sensory organs (Campuzano and Modolell, 1992). In some places, a fixed number of large 
sensory bristles develop at invariant locations, and this was traced to the activity of specific 
cis-regulatory enhancers in the Achate-Scute Complex locus that integrate positional cues 
and direct proneural expression in small clusters of cells (GOmez-Skarmeta et al., 2003). 
Here, Notch-mediated cell-cell interactions are understood to simply refine the readout of 
positional cues. As a caveat, this model does not readily account for the emergence of SOPs 
near the center of proneural clusters, where inhibition is expected to be strongest, unless 
one assumes a sharp spatial bias that favors a subgroup of cells (Troost et al., 2015). In broad 
regions of competence, on the other hand, lateral inhibition should produce irregular arrays 
of isolated SOPs. Indeed, some areas of the body are covered with irregular arrangements of 
small bristles. Such "salt-and-pepper" patterns arise in mathematical models of lateral 
inhibition patterning, centered on the dynamics of Notch and Delta, through the 
amplification of small differences between cells (Collier et al., 1996; Heitzler and Simpson, 
1991). Of note, these models, by contrast with models of heterocyst patterning in 
cyanobacteria, do not require cell-autonomous positive feedback, since signaling is 
restricted to a cell's neighbors and there is no self-inhibition. In other areas of the body, 
however, extended arrays of bristles exhibit regular arrangements. In the medial thorax, for 
instance, five rows of bristles develop on each side of the midline (see the SOP pattern in 
Figure 1D). Each of these rows emerges from a proneural stripe, that could reflect an 
underlying prepattern. It appears however that these stripes are not individually specified by 
specific enhancers of the achaete and scute genes (Corson et al., 2017). Instead, they 
develop in a defined temporal sequence that is suggestive of a self-organized process. 
Supporting this view, analysis of proneural and Notch activity over the course of patterning 
revealed that this sequence is governed by Notch signaling. Before proneural genes are 
expressed, a bimodal gradient of Delta gives rise to two bands of Notch activity in regions of 
intermediate Delta levels, forming a negative template for the expression of Achaete and 
Scute in stripes 1, 3, and 5 (Corson et al., 2017). Concomitant with the onset of Achaete and 
Scute expression in these first three stripes, Delta is expressed in stripe 3 and Notch activity 
becomes dynamic. This expression of Delta does not depend on Achaete and Scute (Parks et 
al., 1997), suggesting that Notch independently serves as a negative template for the 
expression of Delta. As the first three stripes begin to resolve, two bands of low Notch 
activity appear between them, allowing the emergence of stripes 2 and 4. Within each 
stripe, proneural activity is initially higher in the center that on the sides, and SOPs are 
eventually selected from the more central cells, producing rows of regularly spaced SOPs. 


These observations suggested that in the medial thorax, Notch signaling governs both 
proneural patterning and SOP selection - a suggestion that was supported by a simple model 
for Notch-mediated patterning (Figure 1E,F). This minimal model, with a single scalar 
variable for the state of each cell, describes how cell states vary according to inhibitory 
signaling activity; signal production is a function of cell state and propagation is 
instantaneous. Simulated in a two-dimensional field of cells, with an initial pattern of 
signaling activity that mimics the early Delta gradient, the model recapitulated the 
sequential emergence of stripes and their resolution (Corson et al., 2017). The model further 
suggested that cell-intrinsic positive feedback, although not strictly required for lateral 
inhibition patterning, accounts for the emergence of SOPs at the center of the stripes. Such 
feedback, which could result from proneural self-activation and/or Notch cis-inhibition, and 
which gives rise to a bistable response, allows cells with higher proneural activity to evade 
inhibition, even though they are exposed to higher Delta levels. The model also required a 
nonlinear progression of signal production as a function of cell state, which could result from 
the action of signaling modulators such as Neuralized, and shapes the transition from mutual 
inhibition among proneural cells to lateral inhibition from SOPs (Corson et al., 2017). With 
these ingredients, the model encapsulates the multi-step process of proneural cluster 
resolution that was proposed earlier (Troost et al., 2015), with no need for specific 
mechanisms to restrict competence to a subgroup of cells or pre-select SOPs; instead, the 
interplay between cell-intrinsic dynamics and signaling drives a gradual refinement of cell 
fates. 


Taken together, model and experiments indicated that positional cues (a gradient of Delta 
that is decoded through cis-inhibition and trans-activation) and cell-cell interactions combine 
to create a regular pattern of proneural stripes and bristle rows (Corson et al., 2017). In 
other tissues, different spatial or temporal cues impinging on lateral inhibition patterning 
may give rise to different cell fate arrangements. In the Drosophila eye, for instance, 
patterning along a traveling differentiation front produces a regular array of R8 
photoreceptors, arranged in a hexagonal lattice (Lubensky et al., 2011; Roignant and 
Treisman, 2009). While the overall logic of bristle patterning is clear, some issues remain 
open. For instance, how inhibitory signaling extends beyond immediate neighbors, as 
required in the model to account for the observed SOP spacing, remains to be established: 
although long-range signaling may involve basal filopodia (Cohen et al., 2010; De Joussineau 
et al., 2003), direct experimental evidence for Delta signaling at basal filopodia is still lacking. 
Also, while Notch cis-inhibition (del Alamo et al., 2011) has been proposed to yield a sharper 
segregation of cell fates (Barad et al., 2010; Sprinzak et al., 2011), whether cis-inhibition 
actually contributes to SOP selection remains to be tested. In summary, a self-organized 
process guided by temporal and spatial, or boundary, cues orchestrates the patterning of 
sensory organs in the dorsal thorax of Drosophila (Corson et al., 2017). In this context, the 
number of bristles per row varies between animals, and with the size of the field. By 
contrast, variability in the number of rows is much more limited, e.g. starved flies can 


partially miss a late-appearing row whereas well-fed over-grown flies can exhibit an extra 
row (Simpson et al., 1999). In other words, this self-organized process results in a 
stereotyped outcome that is relatively insensitive to physiological perturbations and genetic 
variation. Self-organization guided by stereotyped initial/boundary conditions may thus be a 
path to reproducible outcomes during development. 


A Turing system for digit patterning in tetrapods 

While discrete spacing patterns arising through cell-intrinsic positive feedback and lateral 
inhibition can be conceived as a limiting case of Turing models, these are more commonly 
associated with the specification of larger-scale elements. In particular, reaction-diffusion 
models have been proposed to underlie the development of various body structures in 
vertebrates, including digits in tetrapods (Cooper, 2015) and ridges on the palate of 
mammals (Economou et al., 2012; Lan et al., 2015). The idea that reaction-diffusion may 
account for the patterning of skeletal elements in developing limbs dates back to the 70s 
(Cooper, 2015; Newman and Frisch, 1979; Wilby and Ede, 1975). Since then, experiments in 
chick and mouse have provided a detailed view of the molecular and cellular processes that 
direct the digit/interdigit decision, which now allow theory to be more explicitly confronted 
with experiments (Miura and Shiota, 2000; Onimaru et al., 2016; Raspopovic et al., 2014; 
Sheth et al., 2012). 


Vertebrate limbs appear as small buds produced by the epithelial-to-mesenchymal transition 
of trunk epithelial cells at specific locations along the body (Gros and Tabin, 2014). These 
buds produce a small number of digits, e.g. three in the chick forelimb and five in mouse 
limbs, that are morphologically distinct along the anterior-posterior (a-p) axis. Classical 
grafting experiments indicated that positional information, in the form of a Shh gradient, 
instructs the morphological identity of digits (but see also Delgado and Torres, (2016) for a 
review of the experimental evidence challenging this view). While in principle positional 
information could be interpreted to generate a periodic pattern, mutations in components 
of the Shh signaling pathway result in polydactyly due to an increase in cell proliferation and 
a larger limb (Towers et al., 2008; Zhu et al., 2008) rather than in altered periodicity of the 
pattern. This indicates that Shh regulates the size of the digit field but is not required for 
digit patterning per se. Moreover, loss of positional identity, through dissociation, re- 
aggregation and grafting of limb cells, does not block digit formation (Zwilling, 1964). Thus, 
generating digits and specifying morphological identities are largely independent processes. 


The idea that digit patterning is self-organized, and might be described by a reaction- 
diffusion model, is supported by several observations. First, time-course analysis of the 
emergence of digit condensations in the developing limb indicates that digits appear as 
regularly spaced condensations in a defined temporal order as the digit field grows (Hiscock 
et al., 2017; Raspopovic et al., 2014; Zhu et al., 2008). Second, an allelic series of mutations 


in the Hoxa13, Hoxd11-13 and Gli3 genes produces mouse mutant embryos with a gradual 
increase in digit number, ranging from five up to thirteen. While this phenotype may reflect 
the changing wavelength of a Turing-like pattern of stripes (Sheth et al., 2012), the 
observation that digits may be initiated as spots and extend with growth, rather than being 
specified de novo as stripes (Hiscock et al., 2017), suggested a different interpretation of 
these phenotypes: spots develop with the same spacing along the a-p axis but at a more 
distal position within the limb, where the digit field is broader (Hiscock et al., 2017). In this 
view, digits appear at regular intervals and their number is defined by the size of the digit 
field, as a Turing-like system would predict (Hiscock et al., 2017). Third, using Sox9 as an 
early digit marker, dissociated limb cells could form periodic patterns ex vivo (Raspopovic et 
al., 2014) [see also (Miura and Shiota, 2000; Ros et al., 1994)]. Together, these data support 
the view that a Turing-like system may be at play. 


What then is the molecular dynamics underlying this patterning system? A transcription 
factor, Sox9, and two signals, BMP and Wnt, were recently proposed to be part of a three- 
node Turing system (Figure 2)(Raspopovic et al., 2014). Indeed, Sox9 expression, marking the 
future digits, is in phase with BMP activity (phosphorylated Smad) but out of phase with Wnt 
signaling activity (b-catenin). Loss and/or inhibition of Sox9, BMP and Wnt activities disrupt 
patterning with BMP activating and Wnt repressing Sox9 expression (H. Akiyama et al., 2002; 
Raspopovic et al., 2014). Absent a full characterization of the interactions between Sox9, 
BMP, and Wnt, an exploration of possible network structures identified a minimal topology 
that could support digit patterning. In a fixed domain and with uniform parameters, the 
model produces labyrinthine patterns, but with suitable conditions, incorporating growth of 
the limb and plausible proximo-distal gradients of its parameters, it recapitulates the 
sequence of digit formation, as seen in patterns of Sox9 expression. Consistent with the 
model, combined inhibition of BMP (digit loss) and Wnt (interdigit loss) resulted in fewer and 
larger digits (Raspopovic et al., 2014). In addition, the expression and function of Sox9, BMP 
and Wnt appear to be conserved for the skeletal patterning in the fins of fishes (Onimaru et 
al., 2016). Nevertheless, some questions remain unanswered. First, since none of the Wnt 
ligands is known to be expressed in a periodic manner, it is not clear how Wnt signaling is 
spatially restricted. Second, BMP activity (phosphorylated Smad) and Bmp2 expression, 
which are subsumed under the same variable in the three-node model, are actually out of 
phase. While this observation could be explained by a more complex model, in which Wnt 
downregulates phosphorylated Smad downstream of BMP (Marcon et al., 2016), whether 
and how this interaction occurs remains to ascertained. 


Taken together, the published evidence makes a compelling case that digit patterning 
involves a Turing-like system, yet the precise nature of the relevant interactions, and their 
molecular basis, remains to be clarified. Also, more stringent tests of the model would be 
desirable. Indeed, while the proposed model matches the data, the latter do not necessarily 
exclude other models, e.g. cell-based models such as cell mobility towards cartilage-like 
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condensations (discussed in (Hiscock and Megason, 2015)). If digits appear as spot-like 
condensations, the temporal order of emergence of these spots should be revealing. Thus, 
further tests of the simple three-node Turing model (Raspopovic et al., 2014) and of its more 
realistic five-node version (Marcon et al., 2016) may require an analysis of the dynamics of 
spot emergence upon experimental perturbations predicted to affect the number, size 
and/or spacing of these spots. In the future, combining endogenous GFP-tagged signaling 
molecules and activity reporters with ex vivo live imaging should permit to confront the 
model with experimental data in a quantitative manner. 


Conditions for patterning in mathematical models 

The example of digit patterning, like other systems involving multiple signals and/or 
transcription factors, raises the question of whether self-organized patterning is subject to 
any general principles. The study of diverse patterning models, in particular reaction- 
diffusion models, has suggested that spontaneous pattern formation depends ona 
combination of local activation and long-range inhibition. Patterned cell aggregation can also 
be interpreted in this way, with local enhancement occurring through cell attraction towards 
regions of higher density, and long-range inhibition through cell depletion and/or 
mechanical forces (see below). In the classic activator-inhibitor model, patterning requires a 
greater diffusivity of the inhibitor (Box 1), and in a model where a non-diffusing species 
induces a diffusing activator and a diffusing inhibitor, the diffusion range of the inhibitor 
must be larger (Box 2). But whether local activation and long-range inhibition are always 
required, in Turing or other models, remains unclear, and has been called into question by a 
recent study of patterning in reaction-diffusion models. Marcon et al. (2016) conducted a 
numerical screen for network topologies that can form Turing patterns. In addition to 
networks where patterning depends on the diffusion parameters, as one might expect, the 
screen also identified networks where patterning is independent of the diffusion rates. All 
that can be said about networks that pattern, Marcon et al. (2016) suggest, is that diffusion 
shifts the balance between stabilizing and destabilizing feedbacks. Still, one could make the 
case that the behavior of some of the networks identified as independent of the diffusion 
rates is compatible with the classic interpretation. In one such example, a non-diffusing 
species induces a diffusing activator, which in turn induces a diffusing inhibitor (see Box 2). 
Here, the effective range of inhibition compounds the diffusion of activator and inhibitor, 
thereby exceeding the range of activation regardless of diffusion rates. A more systematic 
investigation may clarify whether the patterning behavior of arbitrary networks can be 
understood along the same lines, or whether some follow an inherently different logic. 


As a possible entry into this question, generalizing on the previous example, one can 
introduce an "interaction function" that describes how the state of a cell varies according to 
the states of surrounding cells (see Box 2). The above-mentioned model for bristle 
patterning in flies (Corson et al., 2017), which treats signaling as instantaneous, assumes 
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such a structure from the start. For a reaction-diffusion model, an interaction function is 
most naturally defined when a separation can be made between variables corresponding to 
an internal "cell state", and to diffusing signals, as in the digit model of (Raspopovic et al., 
2014). Of note, such a separation is absent in the classic activator-inhibitor model. And more 
broadly, it may not be operative if one allows positive feedback among diffusing species, 
since patterning may then occur without changes in cell state. But when cells interact 
through diffusing substances that are produced in a cell-state-dependent manner, it is 
possible to ascribe a signaling field to each cell, and therefore to define an interaction 
function that couples the cell states. And in that case, it seems reasonable to conjecture that 
patterning should depend on a distance-dependent balance between activation and 
inhibition. Beyond reaction diffusion-models, the notion of interaction function may provide 
a common framework to analyze different interaction modalities, whether chemical and/or 
mechanical, and abstract, to some degree, the dynamics of pattern formation from the 
specifics of interaction networks. 


Integrating mechanics with signaling: feather and hair patterning in amniotes 

To date, most investigations of self-organized patterning have focused on the action of 
molecular signals. Yet models involving cell migration and/or mechanical forces can give rise 
to similar patterns, such as spots in two dimensions, and there are multiple instances where 
they have been proposed as plausible alternatives to chemical models, as in the patterning 
of digits, as well as feathers and hair (Murray et al., 1983). These skin structures, which cover 
the body surface of birds and mammals, develop from follicles that form at regular intervals 
during organogenesis (Figure 3A), and a Turing-like mechanism has long been proposed to 
explain their arrangement (Nagorcka and Mooney, 1985; 1982). Recent studies in chick and 
mice have shed new light into the molecular and cellular underpinnings of this patterning 
process (Glover et al., 2017; Shyer et al., 2017), and on the relative contributions of chemical 
signals and mechanical forces. 


The formation of follicles involves communication between two layers of skin cells that are 
separated by a basement membrane: the dermis, composed of mesenchymal cells, and the 
overlying epidermis, a sheet of epithelial cells. At mid-embryogenesis, cells reorganize to 
form regularly spaced dermal condensates, i.e. multicellular aggregates of mesenchymal 
cells in the dermis, associated with closely packed epidermal cells forming placodes. These 
morphological changes are accompanied by changes in gene expression in both skin layers. 
Experiments in which the epidermis and dermis from different body parts were separated 
and recombined suggested that patterning first occurs in the dermis and that patterning 
information is then relayed to the overlying epidermis (Dhouailly, 1975). 


Two types of models have been considered to explain how a pattern emerges in the dermis, 
and how it is possibly integrated with pattern formation mechanisms operating in the 
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epidermis (Painter et al., 2012). In one, a molecular pattern of gene expression emerges 
first, presumably in the dermis, and secondarily directs changes in tissue organization to 
produce the follicles. Here, patterning has a chemical basis, possibly involving a Turing 
system. In a second model, mobile mesenchymal cells first organize into regularly spaced 
aggregates in the dermis, and this pattern of cell aggregates then induces changes in gene 
expression in both skin layers. 


So far, most studies have addressed the first type of model and focused on the identification 
of the relevant diffusive signals operating as activator/inhibitor pairs, such as Wnt/Dkk in the 
mouse (Sick et al., 2006) and BMP7/BMP2 in the chick (Michon et al., 2008). There is reason 
to doubt, however, that any pair of signals can account for skin patterning. Indeed, Shh and 
FGF are also implicated as activators in the chick, and a recent study that examined the 
interactions between BMP, FGF, and Wnt in the mouse, suggests that positive feedback may 
occur only indirectly, through multiple inhibitory interactions (Glover et al., 2017). Thus, the 
precise function of diffusible signals in skin patterning remains unclear. Could it be instead 
that mechanical cues are the main driver of skin patterning? Since Wnt/B-catenin is active 
and required in both the dermis and epidermis for skin patterning in mice (Chen et al., 2012; 
Zhang et al., 2009), and since nuclear B-catenin is a known read-out of mechanical cues 
(Benham-Pyle et al., 2015), it is conceivable that mechanical interactions between 
contractile mesenchymal cells drive the spontaneous formation of cell aggregates in the 
dermis and that these dermal aggregates in turn provide an endogenous source of 
compression for the epidermis. Consistent with this, a recent study in chick indicated that 
nuclear accumulation of B-catenin did not precede the formation of dermal condensates and 
epidermal placodes, and the loss of nuclear B-catenin in drug-treated explants did not 
prevent the formation of multicellular aggregates (Shyer et al., 2017). Also, while untreated 
explants developed regularly spaced follicles, skin explants under compression (enhanced 
actomyosin contractility) showed uniform nuclear B-catenin and broad expression of follicle- 
specific genes. Conversely, nuclear B-catenin and follicle-specific gene expression were lost 
in relaxed explants (under reduced actomyosin contractility). The formation of 
condensations was also sensitive to the stiffness of the substrate on which the explants were 
grown, suggesting that they arise through a mechanical instability that only operates in a 
certain range of mechanical parameters. The precise nature of this instability remains to be 
clarified, likely requiring further experimental studies that directly probe mechanical forces 
within the tissue or compare the effect of substrate stiffness with the application of forces at 
the tissue border, combined with modeling of force generation and transmission in this 
layered system. Nevertheless, these results support a model where mobile mesenchymal 
cells first form a periodic pattern of cell aggregates and this patterning information is then 
relayed to the epidermis through dermis-induced compressive forces (Figure 3B). Thus, in 
the chick, mechanical forces are proposed to drive the patterning of feather follicles, and 
nuclear B-catenin integrates mechanical cues to regulate gene expression in emerging 
follicles (Shyer et al., 2017). 
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Is the same process involved in hair formation in the mouse? A recent study suggests that 
patterning in the epidermis actually precedes the detection of cell aggregates in the dermis, 
implying that the latter may not pattern the epidermis as in the chick (Glover et al., 2017). 
Rather, mesenchymal cells may respond to epidermal patterning cues to form dermal 
condensates in the mouse. Nevertheless, when epidermal patterning is disrupted, as 
monitored by the loss of nuclear f-catenin activity under BMP inhibition and high FGF, a 
regular pattern of dermal condensates is observed. Thus, patterning can take place in the 
dermis in the absence of patterning cues from the epidermis, possibly through self- 
organized aggregation of mesenchymal as in the chick, although it may normally be driven 
by a combination of diffusive signals operating in the epidermis (Figure 3C) (Glover et al., 
2017). Therefore, self-organized cell aggregation (in the dermis) and reaction-diffusion (in 
the epidermis) may converge for hair patterning in the mouse. Interestingly, these two 
mechanisms make different predictions about the positioning of follicles at tissue 
boundaries (Glover et al., 2017). Indeed, any mechanism based on cell aggregation is fueled 
by mobile cells migrating towards the emerging condensates; since the supply of mobile cells 
is limited next to the edge of the tissue, condensates are predicted to form at a distance 
from the edge. By contrast, in reaction-diffusion systems based on local activators and long- 
range inhibitors, the edge serves as a sink for diffusible signaling molecules, hence favoring 
local activation over long-range inhibition so that follicles will tend to form along tissue 
boundaries. These predictions were tested by introducing a cut edge in skin explants and by 
examining the position of dermal condensates relative to the cut edge (Glover et al., 2017). 
In control explants, cell aggregates were found along the edge, consistent with patterning by 
diffusible signals. On the other hand, when explants were treated with BMP inhibitors and 
FGF to inhibit patterning in the epidermis, condensates formed away from the edge, 
consistent with a driving role for cell aggregation. It thus appears that the ability of dermal 
cells to self-organize is conserved in amniotes but, in the mouse, is slaved to patterning in 
the epidermis (Glover et al., 2017). 


Further elaborating on the interactions between dermis and epidermis in skin patterning, a 
recent study in chick has identified FGF expressed in the epidermis as an attractant for 
dermal cell aggregation (Ho et al., 2019). FGF expression was enhanced upon tissue 
compression, suggesting that mechanical feedback on gene expression in the epidermis is an 
integral part of dermal condensate patterning in chick (while mesenchymal Transforming 
Growth Factor B signaling may support autonomous aggregation in mouse (Glover et al., 
2017)). This study also implicated a "priming wave" of Ectodysplasin A (EDA) signaling and a 
tissue-wide gradient of cell density in the sequential formation of the pattern, which is 
initiated at the dorsal midline and progresses outwards. Specifically, experiments in which 
EDA signaling and/or cell density were perturbed, combined with a mathematical model 
introduced by Painter et al. (2018), indicated that EDA modulates a critical cell density above 
which aggregation can occur, such that the spread of EDA activity, together with a time- 
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dependent profile of cell density, control the timing of follicle formation. The seguential 
emergence of follicles is essential to produce a regular arrangement, and this was lost when 
follicle formation was decoupled from the priming wave. As with patterning in the 
Drosophila eye (Lubensky et al., 2011; Roignant and Treisman, 2009), this is an instance 
where local self-organization, combined with a wave that directs sequential patterning, 
result in tissue-wide spatial order. 


The patterning of hair and feathers nicely illustrates how self-driven morphogenesis and self- 
patterning (Sasai, 2013) may combine to produce patterns of cell fates. Indeed, the 
integration of mechanics with chemical signaling is an emerging theme in pattern formation, 
e.g. quorum sensing by Yap that sets a threshold for the oscillatory expression of clock genes 
during somitogenesis (Hubaud et al., 2017); buckling and folding of the gut epithelium (in 
response to compressive stress resulting from smooth muscles resisting the growth of the 
gut) creates pockets of mesenchymal cells where Shh concentrates, creating a regular 
distribution of regions with high Shh (devoid of stem cells, corresponding to the villus tip) 
and regions of low Shh (enriched in stem cells, corresponding to the crypt) (Shyer et al., 
2013; 2015). In this regard, it is interesting to note that several cell-cell signaling pathways, 
e.g. YAP/TAZ, integrin, and Notch, include mechanosensitive molecules. 


Stripe patterning in fishes: a Turing-like system based on interacting cells? 

Most animals exhibit species-specific pigmentation patterns that may serve for camouflage, 
sexual communication, or mimicry. Theory has long been used to account for these patterns 
(Koch and Meinhardt, 1994; Murray, 2003). In a landmark study, Kondo and Asai (Kondo and 
Asai, 1995) proposed that a reaction-diffusion system can account for the stripe pattern 
seen in adult angelfish. Specifically, the pattern changes recorded over a 2-month period of 
body growth, including the intercalation of new stripes, matched computer simulations 
based on a reaction-diffusion model relying on short-range activation and long-range 
inhibition (although the model was simulated in 1D, and may therefore not support the 
same sequence in 2D (Hôfer and Maini, 1996)). These observations prompted the search for 
the underlying fast-diffusing inhibitor and slow-diffusing activator, and zebrafish was chosen 
as a model organism to decipher the molecular and developmental basis of stripe formation. 
The ensuing studies have provided further support for self-organized patterning, like the 
observation that, following a local erasure of the stripes, a labyrinthine pattern regenerates 
and neighboring stripes can be displaced, as predicted by a Turing model (Nakamasu et al., 
2009; Yamaguchi et al., 2007). All the while, they have revealed that it involves a much richer 
repertoire of cell-cell interactions than the simple diffusion of signaling molecules. 


The zebrafish owes its name to a series of longitudinal blue and yellow stripes in the trunk 
and in the fins (Figure 4A). These stripes are produced by the multi-layered arrangement of 
three types of cells present in the skin, the yellow xanthophores, the blue or silvery 
iridophores, and the black melanophores. Patterning has been best studied in the trunk 
where the stripe pattern develops in two separate stages. First, a simple pattern of four 
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stripes of black melanophores forms through the migration of neural crest-derived cells 
during embryogenesis (Singh et al., 2016; Singh and Nüsslein-Volhard, 2015). Iridophores 
appear to be associated with melanophores, whereas xanthophores are mostly scattered. 
This larval pattern is essentially complete after six days of development. Second, the adult 
pattern develops over a one-month period starting three weeks after fertilization. It is in 
part contributed by the larval xanthophores, which first de-differentiate and proliferate to 
cover the trunk as a top layer and then re-differentiate to contribute to the light stripes 
while remaining undifferentiated in the dark stripe regions. By contrast, the larval 
iridophores and melanophores do not contribute to the adult pattern. New iridophores and 
melanophores are produced during metamorphosis by neural crest-derived stem cells that 
are located deep inside the body (Mahalwar et al., 2014; McMenamin et al., 2014). These 
iridophores and melanophores (or their immediate precursors) migrate via specific routes 
towards a skin territory that is largely covered by the yellow xanthophores. Once in the skin, 
iridophores continue to proliferate whereas melanophores gradually stop dividing and 
enlarge (Mahalwar et al., 2014; Singh et al., 2016; Singh and Nusslein-Volhard, 2015). 


How then do the yellow (or light) and blue (or dark) stripes emerge? A simplified view, 
synthesizing a large body of elegant lineage and mutant studies, is presented here (Figure 
4B)(Singh and Nüsslein-Volhard, 2015). The first horizontal yellow stripe is organized by 
iridophores that reach the skin along the myoseptum and interact with the overlying 
xanthophores to promote the production of the yellow pigments. Next, melanophores 
preferentially reach the skin along the sides of this first stripe, where they interact with 
iridophores, which disperse, and xanthophores, which adopt a stellate shape and appear 
faint. As a result, two blue (or dark) stripes appear flanking the central yellow stripe. Then, 
as the juvenile fish grow, additional light stripes form through patterned aggregation of 
iridophores as they delaminate from the light stripes, proliferate, disperse, and re-aggregate. 
The horizontal myoseptum thus seems to serve as a positional cue that orients the stripe 
pattern as it defines where iridophores first aggregate (Singh and Nüsslein-Volhard, 2015). 
Consistent with this, mutants lacking the myoseptum exhibit a labyrinthine pattern of stripes 
with no predominant orientation but normal width (Frohnhôfer et al., 2013). On the other 
hand, what defines the width of yellow stripes and what triggers iridophores to condense to 
form the flanking yellow stripes, hence defining the width of the black stripes, is not well 
understood. Cell-cell interactions are clearly important as revealed by mutant analysis 
(reviewed in (Singh and Nüsslein-Volhard, 2015; Watanabe and Kondo, 2015a)). Mutants 
lacking any one of the three types of pigment cells show strong patterning defects, 
indicating that no single type is indispensable for patterning per se and that all three cell 
types contribute to patterning. The analysis of double mutants, in which only one cell type 
remains, indicates that each type of pigment cell has the ability to cover the entire field 
(Frohnhôfer et al., 2013). This excludes the possibility that pigment cells interpret a fixed 
pre-pattern and instead suggests that patterning involves mutual inhibitory cell-cell 
interactions. The nature of these interactions has been studied in vivo, using laser ablation 
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(Nakamasu et al., 2009), mutant mosaic analysis (Maderspacher and Nüsslein-Volhard, 
2003), altered gene expression (Patterson et al., 2014) and conditional inactivation (Parichy 
and J. M. Turner, 2003), as well as in vitro (Inaba et al., 2012). Together, these studies 
indicate that iridophores attract xanthophores and promote their differentiation, and that 
iridophores and xanthophores repel melanophores at short range and/or promote their 
death. Melanophores may be further cleared from the yellow stripes through long-distance 
signaling via cellular projections produced by undifferentiated xanthophore cells located 
within the black stripes (Eom et al., 2015; Eom and Parichy, 2017), while projections 
extending from melanophores to xanthophores have been proposed to support 
melanophore survival and govern the width of the black stripes (Hamada et al., 2014) - 
although the latter projections may be more prevalent in the adult (Eom et al., 2015). 
Finally, iridophores promote melanophore survival via a yet unknown mechanism 
(Frohnhofer et al., 2013). In summary, long- and short-range cell-cell interactions impacting 
on cell proliferation (iridophores and xanthophores), cell dispersal (iridophores) and 
survival/differentiation (all cells) appear to form the basis of a self-organized patterning 
process. Positional cues, provided notably by the myoseptum, may impinge on this self- 
organized process by specifying the sites where new cells are inserted, resulting in a 
reproducible outcome. 


These developmental studies therefore clearly establish that stripe patterning does not 
involve a simple reaction-diffusion system with a fast-diffusing inhibitor and slow-diffusing 
activator produced by static and equipotent cells. Actually, any model based ona 
stripe/inter-stripe fate choice taken by immobile and equivalent cells would not be relevant 
since the elements of the pattern, i.e. the three types of pigment cells, are mobile and 
display cell-type specific behaviors. In other words, the reaction-diffusion system that was 
initially proposed (Kondo and Asai, 1995) does not have simple molecular correlates (see 
(Singh et al., 2015) for discussion), and is best interpreted as an effective description of 
interacting cell populations. This was made explicit in a later model, with three variables for 
the densities of melanophores and xantophores, and a signal mediating their interactions 
(Nakamasu et al., 2009). The model retains the mathematical form of a Turing model, with 
diffusion substituting for interactions between distant cells. Such a model can recapitulate, 
qualitatively at least, the sequential emergence of stripes from an initial template (Figure 
4C), and the regeneration of a disordered pattern following the laser ablation of skin cells 
(Nakamasu et al., 2009). And a parallel could be drawn between the amplification of small 
perturbations in the model, resulting from mutual inhibition between the two cell types, and 
the "priority effect" whereby the earliest pigment cells to differentiate at a given location 
influence the makeup of the cell population at later times (Patterson et al., 2014). However, 
because the key role of iridophores in stripe patterning, and the complexity of the relevant 
cellular behaviors (proliferation and dispersal, elimination, survival and differentiation to 
cover a large area, or extend long cytoplasmic processes) remain hidden, the relevance and 
predictive value of Turing-like models have been disputed (Mahalwar et al., 2014; Singh et 
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al., 2015; Watanabe and Kondo, 2015b). Recently, a model that explicitly incorporates the 
three cell types has been proposed (Volkening and Sandstede, 2018). This model involves 
both short- and long-range interactions governing a set of behaviors (death, change in 
morphology) for all three cell types. Simulations of these discrete-time rules on a mixed 
population of cells that divide, migrate and interact in an expanding 2D field were used to 
suggest that iridophores contributes to patterning robustness (Volkening and Sandstede, 
2018). More generally, models that explicitly incorporate the three cell types can readily be 
compared with a broader set of experimental observations, e.g. patterning dynamics in 
mutants with only one or two pigment cell types present (Frohnhôfer et al., 2013), and could 
help clarify the extent to which complex cellular interactions can be reduced to a simpler, 
effective description, at least in wild type. In addition, the picture that emerged from studies 
in zebrafish, where positional cues steer a self-organized system towards a defined outcome, 
could be challenged by studies in other species. Within the genus Danio alone, only zebrafish 
show continuous, longitudinal stripes, which can be traced to the horizontal myoseptum. 
Other species exhibit diverse patterns, e.g. vertical stripes (Singh and Nüsslein-Volhard, 
2015), suggesting that alternative directional cues may be at play. 


Evolving skin pigmentation patterns in lizards: a cellular automaton emerging from 
reaction-diffusion systems 

The systems covered so far emphasize individual cells as the elementary units of tissue 
patterning. A recent study of skin patterning in lizards, on the other hand, suggests an 
instance where patterning is best understood as the product of interactions between 
supracellular units. In ocellated lizards, the skin color pattern seen at birth evolves during 
the life of the animal. In juveniles, regularly spaced white ocelli are observed over a brown 
background, with each ocellus comprising a dozen of scales. This pattern gradually changes 
over a three-year period to form a labyrinth-like pattern with about two thirds of the scales, 
initially white or brown, becoming black and the remaining third becoming green (Figure 
5A). During this process, about a third of the scales individually switch color, from green to 
black or black to green; these color changes are not random but rather tend to create strings 
of green or black scales, producing a pattern resembling a labyrinth (Manukyan et al., 2017). 
This remarkable observation begs two questions. First, how is color determined at the level 
of individual scales? Each scale is produced by a dermal condensation that comprises 
pigment-containing chromatophores and light-interfering iridophores. Although the 
molecular basis of scale color specification is not well understood, it likely involves 
interactions between structural and pigmentary colors (Saenko et al., 2013). In analogy with 
the reaction-diffusion system proposed for stripe patterning in fishes, it is proposed that a 
reaction-diffusion system might regulate the adoption of a single color by a given scale 
through non-linear interactions between chromatophores and iridophores. Second, how 
does a scale switch color as a function of the state of its neighbors? (Manukyan et al., 2017) 
showed that the dynamics of the pattern could be captured by a probabilistic cellular 
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automaton, where each element (scale) admits a discrete state (color) (Figure 5B). At every 
time step, the probability that a scale switches color is a function of its current state and the 
states of its neighbors. Iterated application of this rule, using the experimentally observed 
switching probabilities, produced patterns that were statistically similar to patterns in the 
adult. Tying together the cellular and supracellular spatial scales, Manukyan et al. (2017) 
further proposed that the same cell-cell interactions that lead to a uniform color within 
scales could account for the coupling between scales, if they act - although more weakly - 
across scale boundaries (Figure 5C). Indeed, simulations of a reaction-diffusion model with 
reduced diffusivity at scale borders reproduced the behavior of the cellular automaton made 
of skin scales. Thus, a single self-organized process may operate at two spatial length scales 
that are set by a morphological boundary affecting cell-cell interactions, to form a discrete 
pattern on a supracellular scale. 


Perspectives 

As illustrated in this review, the study of diverse model systems provides unambiguous 
evidence for self-organization in tissue patterning, and the relevant interactions are now 
well enough characterized that quantitative comparisons with mathematical models are 
becoming possible; in addition to interactions mediated by chemical signals, evidence for a 
role of mechanical forces is accumulating, which should stimulate a re-examination of 
mechanical models for patterning. All the while, it also evident that multiple, possibly 
redundant interactions are at play in any given system. And because alternative models 
centered on different interactions (different combinations of signals, or mechanical forces 
vs. Signals) can produce similar patterns, whether we have reached a predictive 
understanding of pattern formation in these systems is unclear. In particular, it could be said 
that a definitive example of a Turing system in development remains elusive. 


How then can we ascertain our understanding of self-organized patterning? Part of the 
answer will come from further quantitative experiments that monitor the dynamics of 
patterning and put models to more stringent tests, such as controlled space- or time- 
dependent perturbations, using e.g. optogenetic and microfluidic tools (Duncombe et al., 
2015; Guglielmi et al., 2016; Tischer and Weiner, 2014). But what we consider a satisfactory 
model also merits consideration. Detailed models can be difficult to define with confidence 
and falsify, and may not provide the most transparent representation of the underlying 
"logic". Maybe one should envision a hierarchy of models, incorporating different levels of 
detail, and embodying different levels of understanding. At the coarsest scale are models 
that describe pattern formation in the form of geometric rules, such as have been proposed 
to account for spacing of hair follicles in mouse (Cheng et al., 2014), or, in the context of self- 
organized morphogenesis, the branching of epithelial structures (Hannezo et al., 2017). Such 
models, while agnostic about mechanism, are readily confronted with quantitative data: 
Where is the next follicle formed, in wild-type or perturbed conditions? How are intervals 
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between heterocysts distributed in Anabeana? And they could be used to constrain finer 
models. Closer to a dynamic picture are models in which interacting elements (typically cells, 
but lizard scales are a supra-cellular example) are endowed with a phenomenological 
"internal state" and coupled through an interaction function (Corson and Siggia, 2017; 
Francois and Siggia, 2012). For instance, one might seek to describe limb or stripe patterning 
with a single state variable, capturing the digit/interdigit or stripe/interstripe decision. In 
these terms, it might be possible to identify conditions on cell-cell interactions that are 
required to recapitulate patterning, independent of how they are realized molecularly. Such 
a description could be challenged by monitoring patterning under aberrant boundary 
conditions, e.g. altered growth of the limb, or following perturbations that are localized in 
space or time. Ultimately, live imaging of gene expression and signaling activity, in both wild- 
type and perturbed conditions, combined with genetic and molecular studies, should make it 
possible to tie more explicitly the dynamics of fate specification to networks of molecular 
interactions. A convenient setting for this may be reconstituted in vitro systems, such as 
have been established to study the self-organization of the early mammalian embryo (Etoc 
et al., 2016). By combining singe-cell RNA sequencing with computational approaches, it is 
also now possible to identify transitional cell states and to uncover developmental trajectory 
of heterogenous cell populations in the embryo (Briggs et al., 2018; Farrell et al., 2018; 
Wagner et al., 2018). These emerging approaches, possibly coupled with lineage tracing 
based on DNA barcodes (Fletcher et al., 2018; Raj et al., 2018; Wagner et al., 2018), will be 
useful to study fate dynamics during patterning, as illustrated by the recent identification of 
dermal condensate precursors during the patterning of the hair follicles in the mouse (Gupta 
et al., 2019; Mok et al., 2019). Finally, experimental progress in the dissection of mechanical 
forces acting within tissues should clarify their function in patterning. 


In a developing embryo, self-organization never operates in isolation. Boundary conditions, 
defined by the geometry of a tissue, extrinsic mechanical forces or pre-existing biases 
(developmental history), may steer patterning towards a defined outcome, or define the 
initial breaking of symmetry. Likely equally relevant in defining boundary conditions are the 
temporal cues that define both the onset and the pace of self-organized processes, but 
these remain largely unexplored, as is the temporal coordination of patterning with other 
developmental processes in the embryo. More broadly, rather than simply upstream of one 
another, positional cues and self-organization may be interlinked (Green and Sharpe, 2015), 
as in bristle patterning where both operate through Delta (Corson et al., 2017). The interplay 
between multiple, overlapping patterning processes may participate in a reliable outcome. 
The precision and robustness of development, once theoretical or speculative 
considerations, can now be approached quantitatively in experiments. Where patterning is 
directed by positional cues, information theory provides a framework to analyze variability 
in gene expression, as in the Drosophila segmentation hierarchy, downstream of maternal 
gradients (Dubuis et al., 2013; Petkova et al., 2019). Whether the same concepts carry over 
to patterns that are refined through cell-cell interactions, or self-organized, remains to be 
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explored. And may require that we embrace the suggestion that information is elaborated, 
rather than simply read out, during development (Oyama, 2000). 
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Box 1 - Spontaneous pattern formation in Turing models 

In a Turing model, static cells interact through diffusing chemicals. Commonly, the diffusion 
ranges of the chemicals are taken to be larger than the scale of individual cells, and their 
interactions are described by partial differential equations for continuously varying 
concentrations. The simplest such model that supports patterning comprises two species. In 
the activator-inhibitor model (sketched below), self-activation of A competes with a 
stabilizing negative feedback through the inhibitor B. Equivalently, in a substrate-depletion 
model, negative feedback occurs through the consumption of a substrate that is required for 
the production of the activator. The equations for an activator-inhibitor model in one 
dimension can be written as 


0,A=f,(A,B)+D,0A 


| 
0 B= f,(A)-k,B+D,0°B k 


where the functions f, and f, describe the interactions between A and B within each cell, 


D, and D, are their diffusion rates, and k, stands for any negative feedback on B, e.g. 


constitutive degradation or self-inhibition. 

One assumes that when the cells are uncoupled (in the absence of diffusion), negative 
feedback dominates, such that the tissue tends towards a stable state with uniform levels of 
the two chemicals. Nevertheless, in the presence of diffusion, it is possible for a pattern to 
spontaneously emerge through the amplification of small perturbations. In a linear stability 
analysis (Cross and Greenside, 2009; Murray, 2003), the time evolution of a perturbation is 
approximated by a linear equation, e.g. for the activator-inhibitor model 


ð, A=@A-fBB+D,0-A 


2 (2) 
ð, B= A-k,B+D,0°B 


where A and B now denote the deviation from the unperturbed state. Diffusion is irrelevant 
(stability is unaffected by coupling) on the largest spatial scales, and dampens small-scale 
perturbations. Perturbations at intermediate length scales, on the other hand, can grow, 
provided that the inhibitor has a greater diffusion length than the activator, where the 
diffusion lengths /, and /, are given by 1. = D,/a@ and IA = D, / kp (Cross and Greenside, 
2009). Intuitively, at length scales that are larger than 7, but smaller than /,, variations in 
inhibitor level are damped, thus self-inhibition is reduced, while self-activation of A is 
preserved. Because stability in the absence of diffusion requires k, > œ, the above condition 
on diffusion lengths implies the condition that is most commonly associated with patterning 
in Turing models, that the diffusivity of the inhibitor is larger than that of the activator. The 
latter condition is more practical since the diffusivities are experimentally accessible physical 


parameters, whereas the “diffusion length" Z, involves the regulation of A. 
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Starting from a uniform background, amplification of perturbations leads to a periodic 
pattern, with a spatial scale that is typically set by the fastest growing perturbation (other 
outcomes are possible, which include traveling waves in models with at least three species, 
but steady waves are most relevant here). In a one-dimensional system, the resulting 
patterns show a regular alternation of high and low activator levels (the figure above shows 


a pattern obtained with the parameter values ~@= B= k,=D,=1and D, =4, such that 


l =land /, =2, and with a stabilizing term — A’ added to the equation for A to produce a 


finite-amplitude pattern). In two dimensions, generic outcomes include stripes and spots. In 
a weakly nonlinear analysis (in the regime close to the instability threshold, where the 
amplitude of the pattern is small), the pattern is described as a superposition of waves of 
different orientations, e.g. the overlap between three stripe patterns at 1202 angles forms a 
hexagonal array of spots. Such an analysis suggests that if the system is symmetric with 
respect to the sign of the perturbation, stripes are selected, whereas an asymmetry favors 
spots (Cross and Greenside, 2009; Ermentrout, 1991). Intuitively, in the symmetric case 
there are equal areas of high and low activator, and they will be arranged in stripes. If we 
consider instead an asymmetric case, e.g. a faster than linear increase in inhibitor vs. 
activator level, the balance between activation and inhibition shifts as the perturbation 
grows, yielding a smaller activated region, arranged in spots. 


Of note, Turing models as described above involve several idealizations. For instance, tissue 
that is truly three-dimensional is commonly described in two-dimensions. While the 
geometry of the tissues of interest, such as epithelia or skin, make this a reasonable 
simplification, it is also possible that exchange of signaling molecules with the surrounding 
medium may contribute to their apparent degradation rate, and/or to their transport. In 
some cases, like angelfish stripes (Kondo and Asai, 1995), the geometry of the pattern may 
motivate a one-dimensional model, yet the stability of the same pattern in a two- 
dimensional setting should be considered (Höfer and Maini, 1996). In addition, what appears 
as simple diffusion in the models reflects complex transport processes, that may be affected 
by the tortuosity of the intercellular space, and hindered or facilitated through binding with 
immobilized molecules or shuttles, and whose regulation in space and time may play 
important roles in patterning (T. Akiyama and Gibson, 2015; Müller et al., 2013). 


Discrete models 

Models that explicitly incorporate discrete cells are required in some contexts, such as the 
emergence of salt-and-pepper patterns through lateral inhibition. In a discrete model, 
ordinary differential equations describe the evolution of concentrations in individual cells. 
Generally, these models lend themselves to the same mathematical treatment, e.g. linear 
stability analysis, as continuous models (Plahte, 2001). All the while, they exhibit some 
specific features: whereas some form of short-range positive feedback, or diffusion, is 
required to maintain smooth concentration profiles in continuous models, cell-intrinsic 
positive feedback (e.g. gene self-activation) and a single, inhibitory signal can support the 
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emergence of a discrete pattern (below left). This can be likened to an activator-inhibitor 
model where the scale of the pattern tends to the scale of a cell, and could describe the 
early stages of heterocyst patterning in Anabaena (cf. Figure 1). In the specific case of 
paracrine signaling, as in Notch-mediated lateral inhibition, an inhibitor alone is sufficient for 
patterning. In the absence of self-inhibition, cell-intrinsic positive feedback is not required. 
Instead, mutual inhibition between neighboring cells acts as a double-negative feedback 
loop, amplifying differences between cells (below right) (Collier et al., 1996). 


Box 2 - The notion of interaction function 

The activator-inhibitor model, with just two species representing both the states of cells and 
their interactions, is an obvious idealization. The simplicity of the model, and its relatively 
intuitive behavior, balancing local activation and long-range inhibition, explain its appeal. Yet 
most patterning processes involve multiple signaling pathways, and an interplay between 
diffusing molecules, such as extracellular ligands, and cell-bound species, such as cell surface 
receptors and transcription factors. Models that explicitly account for this complexity are 
less straightforward to analyze systematically, and whether they lend themselves to an 
intuitive interpretation is unclear. The notion of "interaction function" (Hiscock and 
Megason, 2015) may provide an element of answer. 


As a simple example, consider a model with one immobile species, e.g. a transcription factor, 
C, and a diffusible activator/inhibitor pair, A, B (see sketch below; dashed circles denote 
diffusing species), described by the following equations: 


0, A= f (C)-k,A+ D, 9A 
d,B= f,(C)-k,B+D,0°B (3) 
0, C= f,(A,B,C) 

This is logically equivalent to the model of (Raspopovic et al., 2014) for digit patterning, if 


one associates high C with the interdigit fate (low Sox9; cf. Figure 2). The linearized 
equations for a small perturbation are of the form 


ð, A=C-k,A+D,0°A 
ð, B=C-k,B+D,0.B (4) 
0.C=aA-—BB+yC 


where, as in Box 1, we use the same notation, A, B, C, to denote the perturbation, for 
simplicity. 


The analysis is greatly simplified, and the patterning behavior of the model becomes more 
transparent, if one considers a regime where the diffusing signals are close to a steady-state 
profile, with production, diffusion, and degradation balanced. This approximation is justified 
when the system is marginally unstable and the pattern develops slowly. In that limit, 
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signaling is effectively instantaneous, and can be subsumed under an interaction function 
that couples the dynamics of the immobile species C at different positions (in different cells). 
For a one-dimensional system or perturbation, a cell at position X gives rise to a signaling 
field with profiles of A and B that decay exponentially with the distance x-X 


p sea 
K (x-X)= e | 
NA DK. ; 
| est (5) 
K,(x- X)= 7 
ge OE oe 


per unit of C in the cell (see plots below); Z, and /, are the diffusion ranges of the activator 
rs De ME 

and inhibitor, [= D,/k, and 1, = D,/k,. 

Summing over all the cells (integrating over X) and inserting into the equation for C 

(equation 4) allows the dynamics to be expressed in terms of the immobile species alone, 


0,C(x) = yC(x)+ | K- X)C(X) dX (6) 


In this equation, the influence of a cell at position X on a cell at position x it captured by the 

interaction function (plotted below) 

S |x-X| Jaa] 

K(x- X)=a@K ,(x— X)+ BK, (x- X)= ela e ” 7 
( ) aC )+ BK, ( ) JET >k I (7) 


AA BB 


Equations of this form, which could equally describe signaling mediated by cellular 
projections, have also been considered to describe patterns of activity in neural fields 
(Amari, 1977). Analysis of equation 6 shows that patterning requires /, >/,, i.e. a greater 
diffusion range of the inhibitor. Of note, stability towards small-scale perturbations, and a 
smooth pattern, require that cell-intrinsic dynamics be stabilizing, i.e. Y <0; the digit model 
of (Raspopovic et al., 2014) corresponds to the marginal case Y = 0, and yields a singular 
Sox9 profile, with an infinite slope at the digit-interdigit boundary (see Figure 2C). 

The right-hand plot below shows a pattern obtained with the parameter values 
a=B=1/4, y=-05,k,=k,=D,=1,and D, =4, such that 1, =1 and /, =2, and witha 


stabilizing term —C* added to equation for C; the magenta curve shows the interaction term 
a@A-— BB appearing in that equation. 


Interaction functions can also be defined for motile cells (Hiscock and Megason, 2015), with 
model variables now including cell densities. In addition, different networks of interactions 
can give rise to qualitatively similar interaction functions. As a specific example, the 
interaction function for the network represented below, from (Marcon et al., 2016), has the 
same form as for the network shown above, i.e. 
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K(x)= e 4-—___ 2 8 
) 2k 1, 2k,k, LI S 


The same profile, with short-range positive feedback and long-range negative feedback, can 
be recovered with the diffusion ranges of the activator and inhibitor being swapped, as 
illustrated below (here, @= B=3/4, D,=4 and D, =1, such that /, =2 and J, =1; other 


parameter values are the same as above; notice that A and B have different profiles, but the 
same interaction function and the same interaction term @A — BB are recovered). 


As an implication, experiments that probe the effective interaction between cells may not be 
sufficient to discriminate between candidate networks. On the other hand, if enough 
information is obtained on the structure of the interaction function (e.g. by relating gene 
expression and signaling activity patterns), is may be possible to describe the dynamics of 
patterning without a comprehensive characterization of molecular interactions. Taken 
together, interaction functions may provide a unifying framework to analyze seemingly 
disparate patterning models. 
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Legend of Figures 


Figure 1. Lateral inhibition patterning in one and two dimensions (A-C) Heterocyst patterning 
in nitrogen-deprived cyanobacteria. (A) Expression of HetR-GFP in differentiating heterocysts 
in Anabaena, 7 hours after nitrogen step-down [from (Corrales-Guerrero et al., 2015)]. 
Notice the elevated HetR-GFP expression in groups of adjacent cells. (B) Heterocyst 
patterning involves the self-activating transcription factor HetR, and lateral inhibition 
mediated by two short diffusible peptides, PatS and HetN. (C) Following nitrogen 
deprivation, HetR-expressing heterocysts (green) are selected through PatS-mediated lateral 
inhibition (blue). Inhibition around differentiated heterocysts is later sustained by HetN 
(red). During subsequent growth of the filament, new heterocysts intercalate between 
existing ones, maintaining a regular spacing. (D-F) Sensory organ patterning in the 
Drosophila thorax. (D) In the pupal thorax, sensory organ precursors (SOPs; labeled by 
expression of the transcription factor Sensless) that give rise to bristles in the adult fly are 
arranged in regular rows. (E-F) A simple mathematical model, with a single variable per cell 
standing for proneural activity (red), recapitulates SOP patterning. The model incorporates 
self-activation and Notch-mediated lateral inhibition (green). With boundary conditions that 
mimic early Notch activity in the thorax, the model produces a sequence of proneural stripes 
that resolve into SOP rows, as observed [adapted from (Corson et al., 2017)]. 


Figure 2. A Turing model for digit patterning in vertebrates. (A) Expression of Sox9-EGFP 
marks digits in limb culture [from (Raspopovic et al., 2014)]. (B-C) A proposed Turing 
network for limb patterning, comprised of the transcription factor Sox9 and two diffusing 
species, Wnt and Bmp. Simulated in one dimension, the model gives rise to patterns in 
which regions of high Sox9 alternate with regions of high Wnt and Bmp [after (Raspopovic et 
al., 2014)]. 


Figure 3. Coupling between gene expression patterning and mechanics in feather and hair 
patterning. (A) The emergence of feather and hair primordia [here, feather buds in a chick 
embryo from (Painter et al., 2012)] is marked by changes in gene expression and tissue 
morphology. (B) In chick, it proposed that patterned cell aggregation in the dermis (green), 
through compression of the epidermis (red), drives gene expression patterning (Shyer et al., 
2017). (C) In mouse, while dermal cells retain the ability to self-organize, this may be 
normally preceded and driven by gene expression patterning in the epidermis (Glover et al., 
2017). 


Figure 4. Stripe patterning through interactions between pigment-carrying cells in zebrafish. 


(A) The stripe pattern in an adult zebrafish (image courtesy of Laure Bally-Cuif). (B) 
Schematic of cell-cell interactions involved in stripe patterning. The first light stripe forms 
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through the interaction between iridophores reaching the skin atthe myoseptum and 
xanthophores. Interactions between xanthophores and melanophores, which are inhibitory 
at short range and activating at long range (dotted arrow), give rise to the adjacent dark 
stripes [after (Singh and Nüsslein-Volhard, 2015)]. (C) A Turing-like model based on the 
interactions between xanthophores (yellow) and melanophores (black), with an initial 
template comprised of one central stripe, produces a sequence of stripes that progresses 
outwards from the midline [after (Nakamasu et al., 2009)]. 


Figure 5. A multiscale model for skin color patterning in lizards. (A) Adult ocellated lizards 
exhibit a labyrinthine skin pattern, comprised of green and black scales. (B) A cellular 
automaton model, in which individual scales switch color according to the states of their 
neighbors, produces statistically similar patterns. (C) Interactions between pigmented- 
carrying cells could account for both the adoption of a single color within each scale, and, 
acting more weakly through scale borders (dotted arrows), the interactions between scales 
that drive tissue-wide patterning. (A) and (B) from (Manukyan et al., 2017). 
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